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Abstract 

We present new precision measurements of the y(2S) total and partial widths from exci- 
tation curves obtained in antiproton-proton annihilations by Fermilab experiment E835 at 
the Antiproton Accumulator in the year 2000. A new technique of complementary scans 
was developed to study narrow resonances with stochastically cooled antiproton beams. 
The technique relies on precise revolution-frequency and orbit-length measurements, while 
making the analysis of the excitation curve almost independent of machine lattice parame- 
ters. We study the y{2S) meson through the processes pp —> e + e~ and pp — > J/y + X — > 
e + e~ +X. We measure the width to be T = 290±25(sta) ±4(sys) keV and the combi- 
nation of partial widths T e + e -Tp p /T = 579 ± 38(sta) ± 36(sys) meV, which represent the 
most precise measurements to date. 
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1 Introduction 

A precise measurement of the excitation curve of narrow charmonium resonances 
depends both on the detection technique (event statistics, detector efficiency) and 
on the properties of the beam-energy spectrum. In e + e~ annihilations, the beam- 
energy spread is substantially larger than the resonance width. The BES Collabora- 
tion at BEPC published two measurements of the \y(2S) width in e + e~ collisions at 
center-of-mass energies between 3.67 GeV and 3.71 GeV [1,2]. The combination 
T e + e -Tp p /T was recently measured by the BABAR Collaboration at PEP-II using 
initial-state radiation [3]. In pp annihilations, the event statistics are lower, but one 
can take advantage of stochastically cooled antiproton beams, with FWHM energy 
spreads of 0.4-0.5 MeV in the center-of-mass frame, to measure the width directly 
from the excitation curve generated by scanning the beam across the resonance. 
Fermilab experiment E760 measured the widths of the J/\\f and \j/(2S) mesons [4]. 
The uncertainty was dominated by event statistics and statistical fluctuations in 
the beam position measurements. A sizeable systematic uncertainty was due to the 
measurement of the beam-energy spectrum. In this paper, we present results ob- 
tained in pp annihilations by Fermilab experiment E835 from data collected during 
the year 2000 run. A new scanning technique, together with higher event statistics, 
improvements in the beam position measurement and momentum- spread analysis, 
allow us to reach the highest precision to date. 



2 Experimental technique 

The main features of the experiment are summarized here. A full description can 
be found in Ref. [5]. 

In experiment E835, antiprotons circulating in the Antiproton Accumulator inter- 
sect an internal hydrogen gas-jet target. The beam is cooled and decelerated to scan 
charmonium resonances. The operation of the Accumulator for E835 is described 
in Ref. [6]. 
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The E835 detector is a nonmagnetic spectrometer designed to extract, from a large 
hadronic background, electron-positron pairs of high invariant mass as a signature 
of charmonium formation. The apparatus has full acceptance in azimuth, with a 
cylindrical central system and a planar forward system. The central detector system 
includes three segmented hodoscopes, straw-tube and scintillating-fiber trackers, a 
threshold gas Cherenkov counter, and an electromagnetic lead-glass calorimeter. 
In the forward direction, a veto counter and a planar lead-glass electromagnetic 
calorimeter are placed. 

The event selection is described in detail in our paper on \j/(2S) branching ratios [8] . 
The main hardware trigger requires two charged tracks, each defined by a coinci- 
dence between two hodoscope counters, with at least one of the two particles tagged 
as an electron or positron by a signal in the corresponding Cherenkov cell. In ad- 
dition, two energy deposits that are roughly back-to-back in azimuth are required 
in the central calorimeter, with an invariant mass greater than a given fraction of 
the center-of-mass energy. A preliminary off-line selection requires that all e + e~ 
candidates have an invariant mass greater than 2.6 GeV. A maximum-likelihood 
method called 'electron weight' rejects backgrounds, mainly photon conversions 
and Dalitz decays of the pion, that mimic electron or positron tracks in the detector. 
It is based on pulse height in the hodoscopes and in the Cherenkov counter, and on 
the shape of the electromagnetic shower in the central calorimeter. The processes 
pp — > e + e~ and pp — > J/\j/ + X — > e + e~ +X are finally selected using kinematic 
fits requiring a X 2 probability greater than 10~ 4 . The overall efficiency, including 
detector acceptance, hardware trigger, and off-line selections, is about 40% (see 
Section 5), while background contamination is only 0.1% for the e + e~ channel and 
1% for the inclusive channel. 

Two scans of the y/(2S) resonance were performed, in January 2000 (47 hours of 
data taking) and in June 2000 (21 hours). For each run i, the luminosity Jz^ and the 
number of selected events /V,- are shown in Tables 1 and 2. 

The resonance parameters are determined from a maximum-likelihood fit to the 
excitation curve. For each data-taking run (subscript i), we assume that the average 
number of observed events jU; in each channel depends on a Breit-Wigner cross 
section Oewr and on the center-of-mass energy distribution, B t , as follows: 



£i J OBWr ( W) B i (w) dw + O hkg 



(1) 



where w is the center-of-mass energy, £, is the detector efficiency, is the inte- 
grated luminosity, and o^ke is a constant background cross section. The integral 
is extended over the energy acceptance of the machine. The spin-averaged Breit- 
Wigner cross section for a spin-/ resonance of mass M and width T formed in pp 
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annihilations is 



(27+1) 16^ (r in r out /r)-r\ 

ffBw(w) - (25+1)2 w2-4m2r2 + 4(w-M)2' (2) 



m and S are the (anti)proton mass and spin, while Ti n and r ou t are the partial reso- 
nance widths for the entrance (pp, in our case) and exit channels. The Breit-Wigner 
cross section is corrected for initial-state radiation to obtain Oewr [4,7]: 



w/2 b 

<JBWr(w)=b J — f — J o B w(Vw 2 -2kw) (3) 



(w/2)" 

= (2/w) b J dtc m i(^w 2 -2t l l b w), (4) 



where the second form is more suitable for numerical integration and b(w) is the 
semiclassical collinearity factor [7], equal to 0.00753 at the \j/(2S). 

The resonance mass M, width T, 'area' (ri n r out /r) and the background cross sec- 
tion Obkg are left as free parameters in the maximization of the log-likelihood func- 
tion log(A) = Y,i^ogP(/ii,Ni), where P([i,N) are Poisson probabilities of observ- 
ing /V events when the mean is /i. 

Both channels pp — > e + e~ and pp — > J/y + X — > e + e~ +X are fit simultaneously 
to the same mass and width. Each channel is allowed its own area and background 
cross section. 1 



3 Beam energy measurements 



The center-of-mass energy distribution 2?,-(w) is critical for width and area mea- 
surements. We summarize here the concepts that are essential for the following 
discussion. More details can be found in Refs. [4,6]. 

The revolution-frequency distribution of the antiprotons is measured by detecting 
the Schottky noise signal generated by the coasting beam. The signal is sensed by 
a 79-MHz longitudinal Schottky pickup and recorded on a spectrum analyzer. An 
accuracy of 0.05 Hz is achieved on a revolution frequency of 0.63 MHz, over a 
wide dynamic range in intensity (60 dBm). 

1 The 'area' parameter is usually chosen in the parameterization of the resonance shape 
because it is proportional to the total number of events in each channel. It is less correlated 
with the width than the product of branching fractions. 
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The beam is slightly bunched by an rf cavity operating at / cav ~ 1.25 MHz, the 
second harmonic (h = 2) of the revolution frequency. The beam is bunched both for 
stability (ion clearing) and for making the beam position monitors (BPMs) sensitive 
to a portion of the beam. Therefore, recorded orbits refer to particles bunched by 
the rf system, and their revolution frequency is f lf = f cay /h. The bunched-beam 
revolution frequency f lf is usually close to the average revolution frequency of the 
beam. Each orbit consists of 48 horizontal and 42 vertical readings. As a result 
of hardware and software improvements, these readings are much less noisy than 
E760's [4], as discussed later in the uncertainty estimates. 

From the BPM readings and the Accumulator lattice model, we can accurately 
calculate differences AL in the length of one orbit and another. The main system- 
atic uncertainties come from BPM calibrations, from bend-field drifts, and from 
neglecting second-order terms in the orbit length. Using the dispersion function 
from the lattice model, the gains of the high-dispersion BPMs can be measured by 
varying the beam energy at constant magnetic field. They show calibration errors 
between 3% and 15%. Their systematic effect on AL is about 0.03 mm. Contri- 
butions from calibration errors in the low-dispersion BPMs are harder to evaluate, 
but they should be comparable. Bend-field drifts (due to temperature variations, 
for instance) appear in the orbit-length calculation as changes in momentum. For 
the y(2S) scans, their contribution translates into an uncertainty in AL of about 
0.04 mm. The second-order terms depend on the derivatives of the vertical and 
horizontal orbit slope differences with respect to the reference orbit as well as the 
slopes of the reference orbit itself. An explicit assessment of these terms is not pos- 
sible, because there are not enough BPMs to measure slopes everywhere around the 
ring. Under reasonable assumptions, one would get an error of about 0.005 mm. A 
test was performed to estimate the accuracy in AL. The systematic uncertainty is 
evaluated by using a January 2000 \j/(2S) orbit to predict the length of a very dif- 
ferent \j/(2S) orbit of known length from August 1997, when the machine lattice 
was also quite different. The difference between the known length and the predicted 
length is 0.05 mm out of 474 m. Since orbits and lattices for the runs used in this 
analysis are much closer to each other, this is taken as the systematic uncertainty in 
AL from the beam-energy calculation for these runs. 

The absolute length L of an orbit can be calculated from a reference orbit of length 
Lq\ L = Lo + AL. The calibration of Lq is done by scanning a charmonium reso- 
nance (the \y(2S) itself in this analysis) the mass of which is precisely known from 
the resonant-depolarization method in e + e experiments [9]. For particles in the 
bunched portion of the beam (rf bucket), the relativistic parameters /3 rf and are 
calculated from their velocity v rf = f rf ■ L, from which the center-of-mass energy w 
of the pp system is calculated: w lf = w(f Tf ,L) = m^2(l + y rf ). (The superscript rf 
is omitted from orbit lengths because they always refer to particles in the rf bucket.) 
In the charmonium region, this method yields good accuracies on w. For instance, 
dw/df = 113 keV/Hz (38 keV/Hz) and dw/dL = 149 keV/mm (50 keV/mm) at 
the y(2S) (J/yr). 
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For width and area determinations, energy differences are crucial, and they must 
be determined precisely. In our usual data-taking, where we keep the beam near 
the central orbit of the Accumulator, a particular run is chosen as the reference 
(subscript 0). Energy differences between the reference run and other runs in the 
scan (subscript /), for particles in the rf bucket, are simply 

wf _ w $ = W (yrf l + ALt)- w(g, L ) . (5) 



Within the energy range of a resonance scan, these differences are largely indepen- 
dent of the choice of Lq. For this reason, the absolute energy calibration is irrelevant 
for width and area measurements. Only uncertainties coming from AL are consid- 
ered. 

Once the energy wf for particles in the rf bucket is known, the complete energy 
distribution is obtained from the Schottky spectrum using the relation between fre- 
quency differences and momentum differences at constant magnetic field: 

Ap 1 A/ ^ 
P V f 



where rj is the energy-dependent phase-slip factor of the machine, which is one of 
the parameters governing synchrotron oscillations. (The dependence of r\ on beam 
energy is chosen during lattice design, as described in Ref. [6]; the variation of rj 
within a scan can be neglected.) In terms of the center-of-mass energy, 

n w -< ff 



Within a run, rf frequencies, beam-frequency spectra, and BPM readings are up- 
dated every few minutes. Frequency spectra are then translated into center-of-mass 
energy through Eq. 7, weighted by luminosity and summed, to obtain the luminosity- 
weighted normalized energy spectra for each data-taking run. 

The phase-slip factor is usually determined from the synchrotron frequency. In our 
case, this determination has a 10% uncertainty coming from the bolometric rf volt- 
age measurement [5]. At the y{2S), the synchrotron-frequency method yields a 
phase-slip factor r\ = 0.0216 ± 0.0022. 

The resonance width and area are affected by a systematic error due to the uncer- 
tainty in r\ . Usually, the resonance width and area are positively correlated with the 
phase- slip factor. A larger r\ implies a narrower energy spectrum, as described in 
Eq. 7. As a consequence, the fitted resonance will more closely resemble the mea- 
sured excitation curve, yielding a larger resonance width. For our scan at the central 
orbit (stack 1), the 10% uncertainty in r\ translates into a systematic uncertainty of 
about 18% in the width and 2% in the area. 
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4 Complementary scans 



For precision measurements, one needs a better estimate of the phase-slip factor 
or determinations that are independent of 77, or both. In E760, the 'double scan' 
technique was used [4]. It yielded 77 with an uncertainty of 6% at the ty/(2S) and 
width determinations largely independent of the phase- slip factor, but it had the 
disadvantage of being operationally complex. 

Here we describe a new method of 'complementary scans' to achieve a similar 
precision on 77 and arbitrarily small correlations between resonance parameters 
and phase-slip factor; the technique is also operationally simpler. The resonance 
is scanned once on the central orbit, as described above. A second scan is then 
performed at constant magnetic bend field (most of stack 29, runs 5818-5831). 
The energy of the beam is changed by moving the longitudinal stochastic-cooling 
pickups. The beam moves away from the central orbit, and the range of energies is 
limited but appropriate for narrow resonances. 

Since the magnetic field is constant, beam-energy differences can be calculated 
independently of AL, directly from the revolution-frequency spectra and the phase- 
slip factor, according to Eq. 7. A pivot run is chosen (5827 in our case, subscript p). 
The rf frequency of this run is used as a reference to calculate the energy for par- 
ticles in the rf bucket in other runs. These particles have revolution frequency ff 
and the energy is calculated as follows: 

P V h£ ff 



For the scan at constant magnetic field, this relation is used instead of Eq. 5. Once 
the energy for particles at ff is known, the full energy spectrum within each run is 
obtained from Eq. 7, as usual. 2 

Using this alternative energy measurement, the width and area determined from 
scans at constant magnetic field are negatively correlated with 77 . The increasing 
width with increasing 77 is still present, as it is in scans at nearly constant orbit. 
But the dominant effect is that a larger 77 brings the energy points in the excitation 
curve closer to the pivot point, making the width smaller. In the case of stack 29, a 
10% increase in 77 implies a -10% variation in both width and area. 

2 For the constant-field scan, the energy distributions may be obtained directly from the 
pivot energy by calculating w — w^, instead of using Eq. 8 first and then Eq. 7. The two- 
step procedure is chosen because it is faster to rescale the energy spectra than to re-calculate 
them from the frequency spectra when fitting for 77 . Numerically, the difference between 
the two calculations is negligible (less than 0.2 keV). Moreover, the two-step procedure 
exposes how the width depends on 77 . 



7 



The different dependence of the width on r\ in the two separate scans is shown as 
two crossing curves in Figure 1 . (Statistical errors, ±36 keV for both curves, are not 
shown.) The constant-orbit and the constant-field scan can be combined. The result- 
ing width has a dependence on r\ that is intermediate between the two. An appro- 
priate luminosity distribution can make the resulting curve practically horizontal. 
The combined measurement is dominated by the statistical uncertainty (±25 keV, 
in this case; not shown in the plot). 

Moreover, thanks to this complementary behavior, the width, area and phase-slip 
factor can be determined in a maximum-likelihood fit where rj is also a free param- 
eter. Errors and correlations are then obtained directly from the fit. 



5 Results 

Both channels in both scans are fitted simultaneously, leaving the phase-slip fac- 
tor as a free parameter. The energy distributions are rescaled according to Eq. 7 
for the 'constant-orbit' scan and Eqs. 7 and 8 for the 'constant-field' scan. The 
log-likelihood function is log(A) = [log P(fif e ,Nf e ) ±logP(juf ,N? )] . For each 
channel, the mean numbers of events flf e and lif are evaluated according to Eq. 1. 
We monitor the efficiency of each data run and the efficiencies vary less than 0.4% 
over a single scan. Both scans have the same efficiency e ee = 0.413 ±0.015. 
For the J/yz + X channel, the constant- field scan efficiency is £^ = 0.402 ±0.01 1; 
differences in detection efficiency between the two scans are accounted for by the 
parameter (e^/e^). They are due to a different configuration of the tracking sys- 
tem, which does not affect the e + e~ channel. The likelihood maximization was per- 
formed within the R package [10] and crosschecked with the Minuit code [11]. 
The results of the fit are shown in Figure 1, Figure 2, and Table 3. 3 

The fitted value of rj in Table 3 is consistent with the one determined from the 
synchrotron frequency (Section 3). The relative uncertainty in the phase-slip factor 
(6%) is equal to that from the E760 double scans [4]. 

Possible statistical and systematic sources of uncertainty in the width and area 
are considered. As discussed in Section 3, each beam spectrum is a luminosity- 
weighted sum of individual energy distributions within each run. Statistical fluctu- 
ations of the BPM readings produce random variations of the measured AL which 
systematically widen the beam spectrum, making the resonance width narrower. 
The BPM noise is evaluated from portions of runs with no energy drifts and its 
standard deviation is 0.02 mm for both stacks. As a result of hardware and software 
improvements, this is much lower than the E760 value, 0.2 mm [4]. The effect on 

3 The y(2S) mass from Ref. [9] is used for the absolute calibration of Lq. The value of M 
in Table 3 is not an independent measur ement. 
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width and area due to BPM noise is larger for small beam widths and for runs with 
no energy drifts. In the worst case, it translates into a systematic uncertainty of 
< 8 keV in the width and < 2 meV in the area. We do not correct for this system- 
atic, but uncertainties are assigned to the results of 4 keV and 1 meV, respectively. 

The systematic uncertainties in the luminosity (2.5%) and efficiency (3.6%) 
directly affect the area, but not the width. They are added to obtain an uncertainty 
of6.1%or35meV. 

The absolute energy calibration does not influence the resonance width and area. 
Instead, a systematic error in the AL determination has the following effects: it 
shifts all runs in stack 1 ; it shifts stack 29 with respect to stack 1 . The systematic 
uncertainty of 0.05 mm discussed in Section 3 translates into < 1 keV for the width 
and < 1 meV for the area, and it is therefore neglected. No systematic uncertainties 
are therefore introduced by combining the two scans. The pivot run in the constant- 
field scan (Section 4) is taken near the central orbit and has a small AL, so that its 
energy relative to the constant-orbit scan can be calculated accurately from Eq. 5. 

Our final results are the following: 



r = 290±25(sta)±4(sys) keV (9) 
r> e -r> /7 /r = 579±38(sta)±36(sys) meV. (10) 



6 Discussion 



A comparison between this width measurement and those of E760 [4] and BES [ 1 ,2] 
is shown in Figure 3. All three values are consistent. The E835 measurement is 
the most precise. Our measurement of (r e + e -Tp p /r) is also compatible, but much 
more precise, than that reported by BABAR, r e+e -r Pp /r = 0.70±0.17±0.03 eV [3] 

This new method of complementary scans can be applied to future experiments for 
the direct determination of narrow resonance widths in antiproton-proton annihila- 
tions (such as PANDA at the future FAIR facility in Darmstadt). If one performs 
a scan at constant orbit and a scan at constant magnetic field in conditions similar 
to those in the Antiproton Accumulator, the uncertainty is mainly statistical. More- 
over, by appropriately choosing the relative luminosities and energies of the two 
scans, one can make the width almost uncorrelated with the phase-slip factor, as in 
the E835 case discussed in this paper. 
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Run 


Energy 


Luminosity 


e + e events 


J/\j/ + X events 


i 


wf (MeV) 


&i (nb- 1 ) 


N ee 


Nf 


5006 


3687.585 


15.20 





1 


5009 


3687.632 


37.20 


1 


3 


5012 


3687.373 


44.40 





8 


5013 


3687.343 


42.20 





5 


5015 


3687.121 


37.20 





5 


5016 


3687.080 


45.10 


2 


11 


5019 


3686.760 


80.80 


5 


18 


5022 


3686.471 


41.58 


5 


45 


5023 


3686.453 


37.15 


4 


32 


5025 


3686.012 


98.56 


67 


280 


5027 


3685.678 


15.88 


12 


26 


5028 


3685.667 


45.37 


19 


76 


5029 


3685.848 


43.44 


20 


85 


5031 


3685.643 


80.99 


18 


67 


5036 


3685.338 


21.22 


1 


9 


5038 


3685.334 


61.43 


4 


13 






747.72 


158 


684 



Table 1 
Stack 1 data. 
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Run Energy Luminosity e + e events J/y + X events 



i 


vtf (MeV) 


J^-Cnb- 1 ) 


N ee 


Nf 


5818 


3686.674 


77.81 


15 


65 


5819 


3686.701 


50.01 


7 


49 


5821 


3686.422 


79.13 


27 


142 


5822 


3686.427 


40.63 


18 


75 


5824 


3686.126 


78.34 


37 


257 


5825 


3686.138 


55.50 


27 


175 


5827 


3685.922 


78.80 


52 


291 


5828 


3685.922 


68.44 


41 


264 


5830 


3685.633 


79.07 


25 


149 


5831 


3685.643 


48.52 


19 


100 


5833 


3684.455 


79.31 


1 


11 


5834 


3684.450 


78.35 





10 


5837 


3684.451 


78.68 


3 


10 






892.59 


272 


1598 



Table 2 
Stack 29 data. 
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Parameter Value Correlation matrix 

2 3 4 5 6 7 8 

1 M(MeV) 3686.111 ±0.009 0.02 0.35 -0.07 -0.64 -0.23 0.10 -0.17 

2 T(keV) 290 ±25 0.35 -0.29 0.02 -0.07 0.05 -0.52 

3 i> e -r>p/r(meV) 579 ±38 -0.44 -0.48 -0.28 -0.58 -0.27 

4 a bkg (e+0 (pb) 3 ±6 0.09 0.07 0.28 0.17 

5 r/ (10 4 ) 216=1=13 0.54 -0.20 0.17 

6 e*/e£(%) 73 ±4 -0.29 0.02 

7 A A e r a i{;X + -f 5.81 ±0.35 -0.17 

8 o hkg (J/ y + X)(pb) 65 ±19 
[log(A)] max -170.95 
Goodness-of-fit tests: 

log-likelihood ratio 68.8 / 50 (P = 4.0%) 

X 2 I d.o.f. 68.3 / 50 (P = 4.3%) 
Table 3 

Summary of the results of the maximum-likelihood fit. 
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Fig. 1. r dependence on Tj for stacks 1 and 29, and from their combination, when the 
phase-slip factor is a fixed parameter. The result of the global fit with free tj is represented 
by the cross. The value of the phase-slip factor from the synchrotron-frequency measure- 
ment (vertical line) and its uncertainty (gray band) are also shown. 
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Constant-orbit scan (stack 1) 




Constant-field scan (stack 29) 
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Fig. 2. y(2S) resonance scans: the observed cross section for each channel (filled dots); the 
expected cross section from the fit (open diamonds); the 'bare' resonance curves Obw from 
the fit (solid lines). The two bottom plots show the normalized energy distributions Bj. 
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Fig. 3. Recent measurements of the y{2S) width. 
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